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(N . Abstract 



We carry out a perturbative determination of mass dependent renormalization factors and 0{a) 
improvement coefficients for the vector and axial vector currents with a relativistic heavy quark 
action, which we have designed to control rriQa errors by extending the on-shell 0(a) improvement 

o : 

■^J- . program to the case of uiq ^> Aqcd with ijiq the heavy quark mass. We discuss what kind 



of improvement operators are required for the heavy-heavy and the heavydight cases under the 
condition that the Euclidean rotational symmetry is not retained anymore because of the rriQa 
corrections. Our calculation is performed employing the ordinary perturbation theory with the 

■ 

fictitious gluon mass as an infrared regulator. We show that all the improvement coefficients 
• i— i . 

are determined free from infrared divergences. Results of the renormalization factors and the 
improvement coefficients are presented as a function of niQa for various improved gauge actions as 
well as the plaquette action. 
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I. INTRODUCTION 

This paper is the third in a series of publications Q,Q| on a new relativistic approach which 
was recently proposed from the view point of the on-shell 0(a) improvement program. The 
generic quark action, proposed first in Ref. jjj, is given by 



s q = E 



m q(x)q(x) + q(x)^ D q(x) + v^<l(x)^ i D i q(x) 



~q(x)Dlq(x) - -y- £ q(x) D*q(x) 

i 

l -^c E ^q(x)a Qi F iq(x) - ^c B E q(x)(TijFijq(x) 



(1) 



While we are allowed to choose r t = 1, other four parameters u, r s , ce and cb should be 
properly adjusted as functions of rriQa and the gauge coupling constant g, in order to achieve 
the 0(a) improvement for all on-shell matrix elements. In Ref. we determine u, r s , ce and 
cb up to the one-loop level for various improved gauge actions. We now report on the 0(a) 
improvement of the vector and axial vector currents at the one-loop level for the relativistic 
heavy quark action. 

In this paper we first make a general discussion about what kind of improvement operators 
are required from the symmetries allowed on the lattice, in which the Euclidean rotational 
symmetry is violated because of mqa corrections. We consider both the heavy-heavv and 
heavy-light cases, where the light quark is massless for the latter. Following Ref. |2| we 
evaluate one-loop diagrams employing the conventional perturbative method with the use 
of the fictitious gluon mass to regularize the infrared divergence. In the massless case this 
method was successfully applied to the calculation of the renormalization constants and the 
improvement coefficients for the bilinear operators^, 0]. 

This paper is organized as follows. In Sees. II and III we determine the renormalization 
constants and the improvement coefficients for the vector and axial vector currents up to 
one-loop level. The results are presented both for the heavy-heavy and heavy-light cases 
as a function of m,Qa with various improved gauge actions in addition to the ordinary 
plaquette action. In Sec. IV we explain how to implement the mean field improvement for 
the renormalization factors. Our conclusions are summarized in Sec. V. Some preliminary 
results are presented in Ref. p. 

The quark and gluon actions and their Feynman rules are already presented in Sec. II 



of Ref. [2|. We employ the notations introduced there without further notice throughout 
this paper. As for the numerical evaluation of the one-loop diagrams relevant for the vertex 
functions of the vector and axial vector currents, we employ the same method used in the 
perturbative determination of the improvement parameters for the relativistic heavy quark 
action, whose technical details are described in Sec. Ill of Ref. Q] • The physical quantities 
are expressed in lattice units and the lattice spacing a is suppressed unless necessary. We 
take SU(iV c ) gauge group with the gauge coupling constant g. 



II. 0(a) IMPROVEMENT OF THE VECTOR CURRENTS 

We consider the on-shell 0(a) improvement of the vector currents both for the heavy- 
heavy and heavy-light cases. Without Euclidean space-time rotational symmetry, the renor- 
malized operators with the 0(a) improvement is written as 

V ***( X ) = zr [q(xh,Q(x) -g 2 c^ d^{q(x)Q(x)} - 9%^ {q(x)Q(x)} 

-g 2 c^{d t q(x)haM x ) - 9 2 c^q(x)lM^Q^)} + 0(g 4 )} (2) 

where Zy tt and Cy^' ' H ' L ^ depend on the quark masses mg and m q . and d~ are defined 
as = <9 M + dfj, and d~ = — . For the time component of the vector currents we can 
choose Cy = Cy = with the aid of equation of motion. In the case of trq = m q we find 
Cy^ = and Cy^ = Cy^ from the charge conjugation symmetry. Once the both quark masses 
are massless, all the improvement coefficients except Cy vanishes. 

In this section Zy tt and CyJ ' H are determined at the one-loop level as a function of 
runa both for the heavy-heavy and heavy-light cases. We employ the relativistic heavy quark 
JL proposed by Z autlfl both for the heav y md it ^ 



A. Determination of the improvement coefficients for the vector currents 

We consider the general form of the off-shell vertex functions of the vector currents on 
the lattice at the one-loop level: 

A£\p,q,m pl ,m p2 ) = lk F* + lk {^ + tf s F*} + {iF\ + i s F^} lk 



3 



+{Pk + q k ) [G k l +i>G k 2 + 4G , i + 
+(p k -q k ) \H k 1+ $H k 2 +iH k ^ 



+ 0{a 2 



(3) 



and 



„,) = 70^1° + 70^2° + ^3° + HM 

+( Po -go) \H 1 +piH 2 +4H 3 + 



0(a 2 



(4) 



with 



3 3 
a=0 a=0 



(5) 



i=i 



where we assume that the off-shell vertex functions are perturbatively expanded as 
A^p, q,m pl , m p2 ) = 7 M + ^(^Y^iP, q,m pl , m p2 ). 



(6) 



The vertex functions (j3J) and are defined for the process depicted in Fig. [T] The coeffi- 
cients F M , G M , if M are functions of p 2 , q 2 , p ■ q, m v \ and m p2 . 

Sandwiching (J3J) and (JD) by the on-shell quark states u(p) and u(q), which satisfy $u(p) = 
im p iu(p) and u(q)<£ = im p2 u(q), the matrix elements are reduced to 



and 



u{q) A k\p, q, m pi> rn p2 )u(p) 

u(qhMp) {F k + im pl F k + im p2 F k - m pl m p2 F k } 

+u{qhkp'su{p) {F k + im p2 F k + im pl F k ] 

+u(q)4slku{j)) [F k + im pl F k + im p2 F k } 

+ (Pk + qk)u(q)u(p) \G\ + im pl G 2 + im p2 G k - m pl m p2 G k ^ 

+ (pk - qk)u(q)u(p) [H k + im pl H k + im p2 H k - m pl m p2 H k } + 0(c 



u{q)A ( o\p, q, m pl , m p2 )u{p) 

u(q)j u(p) {f° + im pl F% + im p2 F° - m pl m p2 F$} 

+ (po + qo)u(q)u(p) [G? + im pX G\ + im p2 G% - m pl m p2 G ^ 

+ {Po - qo)u{q)u(p) [F° + im pl H 2 + im p2 H® - m pl m p2 H®} + 0(c 



(7) 



(8) 



For convenience we express the coefficients as 

X k = + im pl F* + %m p2 F\ - m pl m p2 F%, (9) 

Y k = G\ + im pl G k 2 + im p2 G% - m pl m p2 G\, (10) 

Z k = H\ + im pl H% + im p2 H$ - m pl m p2 H^, (11) 

R k = F* + im p2 Fj + iuipiFg, (12) 

S k = F 5 fc + im pl F* + im p2 F^ (13) 

and 

X = F° + im pl F° + im p2 F° - m pl m p2 F^ (14) 

Y Q = G\ + im p iG 2 + im p2 G° z - nipim^G ^ (15) 

Z Q = H ° + im pl H% + im p2 Hl - m pl m p2 H^. (16) 

Since the above coefficients contain both the lattice artifacts and the physical contribu- 
tions which remains in the continuum, we have to isolate the lattice artifacts in order to 
determine the improvement coefficients in eq.(J2J). The improvement coefficients are given by 





f \r \latt / \r \cont 

— y^k) — y^k) 5 


(17) 


ic v k 


\latt /-tr \cont 

— y^k) ~ y*k) 1 


(18) 


lc v k 


/ry \latt fry \cont 

— y^k) — y^k) , 


(19) 


-< 


= (S k ) latt , 


(20) 


< 


= (R k ) latt , 


(21) 




1 y \latt / y \cont 

— l A oJ - {^0) , 


(22) 


ic v 


/xr \latt / ~\t~ \cont 


(23) 


ic v 


fry \latt / \cont 

— l z oJ - l^oj , 


(24) 



where the continuum contributions are obtained by employing the naive dimensional regular- 
ization (NDR) with the modified minimal subtraction scheme (MS). We have R k = S k = 
in the continuum from the space-time rotational symmetry. Here it is reminded that in case 
of m p2 = m p i we obtain c Vk = Cy k and c Vk = c V[] = from the charge conjugation symmetry. 
The renormalization factor of the vector currents is obtained by 
Z 



latt 

V " ~ ■ (1 " 9 2 A V ) (25) 



^cont ~~ V ^Q,lattA "°plJ\/ ^ 3 ,lattV " l p2 



with 



Z g°Lt( m p2) = cosh(m^) +r 4 sinh( 



,(0) 



cosh(mp°' ) ) + r t sinh(rap° ] 



(0)v 

p2 , 



(0)> 



m 



Ax 



7m 2 2 ' 



(26) 
(27) 
(28) 



where Aq i9 denote the wave function renormalization factors, which are already given in 
Ref. Although we evaluate -Zy° nt in MS scheme with NDR in this paper, the reader may 
be interested in the value defined in MS scheme with DRED. The conversion between these 
two definitions is easily done by the relation 
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Zy u nt (NDR) = Zy° nt (E)RED) _ l g \ 



(29) 



Employing a set of special momentum assignments p = p* = (po = im P i,Pi — 0) and 
q = q* = (go = irn P 2,qi = 0) or g = g^ = (g = —im p2 ,qi = 0), where subscripts s and 
d represent the scattering and the decay respectively, we extract the relevant coefficients 
Xk, Yk, Zk, Rk, Sk for Vu from the off-shell vertex function (J3J): 

1. 



y 1 
1 1. 



-Tr 



Ajj? (1 + 70)7* 



p=p*,q=q* 



X s k + 2 ? m pl ^ 



+ 2im p2 S s k 



1 

4 

1. 



Tr 



-Tr 



in. 

4 



|_A« (1+7o )_|_A« (1+7o ) 7i7i 

|_ A a> (1 + 7o) _|_ A « (1 + 7o)7i7j 
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-i i^fc 

- P=V*,<i=1*s 

i^k 



Al 1} 7fc (l + 7o ) + 2zm pl ^AS ; 1) (l + 7o ) 7i7 fc 

op. 



p=p*,q=q* 
i^k 



Al 1} (l + 7o ) 7fc + 2im p2 -^-A k 1 \l + 7o ) 7fc7i 
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p=p ,q=<is 

1 i^k 



(30) 
(31) 

(32) 

(33) 

(34) 



->p=p ,q=ii 

where superscripts s and d in X^, Yk, Zk, Rk, Sk represent their momentum assignments. On 
the other hand, we employ (JH) to determine X ,Y , Z for Vq\ 



X c n 



-Tr 



d d 
Aq^o ~ im p i—A^\l + 7o ) 7fc + %m P 2— A^\l + 7o ) 7fc 

c>Pk oqk 



p=p* } q=q* 



X d + im pl (Y d + Z d ) - im p2 (Y d - Z d ) 



-Tr 



A 1} (1 + 7o ) + 2zm p2 ^-A 1) (l + 7o ) 7fc 



X s + im pl (Y s + Z s ) + zm p2 (Y s - Z s ) = -Tr [A^(l + 7o ) 



p=p*,q=q* d 



p=p*,q=qi 



(35) 

(36) 
(37) 



where we have used the fact that F k , G k and H k are functions of p 2 , q 2 and p ■ q, so that 



dpi 
dm 



dpi 
dG 11 



p=p ,q=q 



p=p*,q=q* 



9q 
dm 



dpi 



dG» 



dqi 



p=p ,q=q 



p=p*,q=q* 



p=p*,q=q* 







(38) 
(39) 
(40) 



p=p*,q=q* 
with i = l,2, 3. 

Here we briefly explain how to deal with the infrared divergence in the above coefficients 
at the one-loop level. We basically follow the method employed in Refs. . Suppose the 
vertex function A^q at the one-loop level is written as 

d A k 



A 



(i) 

fe,0 



-I kj0 (k,p,q,m pl ,m p2 , A), 



(41) 



U (2tt) 4 ' 

where A is the fictitious gluon mass introduced to regularize the infrared divergence. We 
extract the infrared divergent term as 

Ik,o(k,p,q,m pl ,m p2 ,X) 

I k ,o(k, p, q, m p i, m p2 , A) - I k ,o(k, p, q, m ph m p2 , A) 

ik,o(k,p,q,m p i,m p2 , A), 



-7T (2tt) 4 ' 

d A k 



+ 



* (2tt) 4 
- d 4 k 



-7T (2tt) 



(42) 



where I k fl(k,p, q, m p i, m p2 , A) should have an analytically integrable expression, whose in- 
frared behavior is the same as I k>0 (k,p, q, m pl , m p2 , A). For I k>0 (k,p, q, m pl , m p2 , A) we employ 



In(k, P, q, m p Y , m p 2 , A) 
6(A 2 - k 2 )C F i lo 



1 



1 



'i{4 + Ji) + m p f\{f + k*) + m£? " ' a k 2 + X 2 ' 

where Cp = (N 2 — 1)/(2N C ) denotes the second Casimir of SU(Af c ) group. A domain of 
integration is restricted to a hypersphere of radius A (< tt) for convenience of an analytical 
integration. 



-Mo 



(43) 
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B. Results for the improvement coefficients of the vector currents 



1. Heavy-heavy case 



We obtain A 7fe , Cy and A 70 , c Vo from eqs.(fT7j). (fTHJl and (|7TT - I23jl choosing m 



pi 



m 



p2i 



where the charge conjugation symmetry demands Cy k 



c" k and c v k = c v = 0- In Fi g s - El 



and El we plot A 7fc ,Cy£ and A 



i 70 ,Cy , respectively, as a function of for the plaquette 
and the Iwasaki gauge actions. The rn$ dependence of Ay M in eq. (J2HJ) is also plotted in 
Fig. HJ The solid lines denote the fitting results of the interpolation: 



A 



7* 



A 



To 



A 



V 



A 



7kl m C°)=o 



+ 



'VI 



m ( ?=0 
pi 



+ 



Ej=l 1 'nZ{ m pl' > }' 
Ef=l ^ni" { m il }* 



i=l u ni \" l pl 5 



-i v di 



E 

l + E? i < ,/, 



+ 



Ef=l ^nf { m pl } 



l + E 



fcZ 



A 



7ol m C°> =0 



o) 
pi 



+ 



A 



™m<°>=0 



+ 



'pi 

Ef=l ^nJ{ m pl' > } t 

pi 



1 + Et^SfK?^ 



(44) 
(45) 
(46) 
(47) 
(48) 
(49) 
(50) 



where we assume c¥r 



at m 



(0) 

Pip 



0. We employ A 



A 70 =0.05169(plaqette), 



0.04802(Iwasaki[8j), 0.04595(DBW2jQ]), c Vk = c Vo =0.01633(plaquette), 0.009728(Iwasaki), 
0.004884(DBW2) and A Vk = A Vo =0.1294(plaquette), 0.06279(Iwasaki), 0.02566(DBW2) 



at m 



(o) 
pi 



0. The values of the parameters v^f and 



k,0 



1, . . . , 5) are summarized in 



Tables HI and ITT1 The relative errors of these interpolations to the data are less than a few % 
over the range < < 10. 
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2. Heavy-light case 



We use eqs. (|T7l-IM)l to determine A 



7fc> C V k 



{+,-,H,L) 



and A 7o , Cy^ \ Assuming that the m p2 a 



corrections are negligible, we evaluate the improvement coefficients, except c v and Cy^ \ 
as a function of m^} with mi "* 



'■ P 2 



yd 



0. In eqs.fSJ), dSZD we find that S s k and Y C 

= 0. Therefore one should extrapolate data at non-zero 
0. We however keep mS^ — 0.0001 in our calculation to determine c v and 



are not determined if we set 

(o) + (o) 
m p { to m v { 

Cy^' ~* since the difference between the value at m p2 



(o) 



rn 



(0) 
p2 



0.0001 and the one extrapolated to 
= is less than 1 %. Figures El and El show the m p °i dependence of A 7fc , Cy~ k '~' H ' L ^ and 
A 70 , Cy o ~\ respectively, for the plaquette and the Iwasaki gauge actions. We also plot Ay M 
in Fig. [3 The interpolation denoted by the solid lines are expressed as 



c 



i-H,L) 
Vfe 



A 



70 



'Vb 



'Vb 



A 



VI) 



A 



7fcU£° 1 ) =0 



+ 



J2j=l V nl{ m plY 

i+e!=i^7{™£V 



™<$=o 



+ 



E 



i=l v n i \ m pl 5 



V 5 /• 

Z—ti=l u ni 



1 + E 5 



Av fc 



i=l u di 



(0) n H 

m pl=° 1 



fc( "' H,i) {Tr4?}* 



i=l \ m pl s 



A 



70 I 



'Vb 



+ 



+ E 

EjLl t, nl{ m pl' > }' 

l + ELi*^}*' 

2jj=l y m l m pl 



Ei=l v ni { m pl } 



l + E 

(o; 
pi 



7 ,o+ 



l + E 



i=l W di 



°-{m p °ly 



A 



E 



oz 

i=l u ni 



Vbl m W=o 



(51) 
(52) 
(53) 



(54) 
(55) 
(56) 
(57) 



with vtf and t>i : ° (i — 1, ... ,5) given in Tables ITTT1 and ITVl Here we use the constraint that 



-Vib 



7i ? 

and c 



vb 



at = 0. The data are well described by these interpolations 



within a few % errors over the range < mil < 10. 



P i 
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III. 0(a) IMPROVEMENT OF THE AXIAL VECTOR CURRENTS 



Let us turn to the axial vector currents. The discussion is in parallel with the case of the 
vector currents. The renormalized operators with the 0(a) improvement is given by 



A^ tt ' R (x) = Z l X lt 



g(»7/i75<20) - 9 2 CA^{q(x)^Q(x)} - g 2 c A d il {q(x)^Q(x)} 
g 2 c L A Ad l(1 (x)H llll5 Q(x) - ^ 2 c^g(x) 7M 757,{^g(x)} + 0(g*)} (58) 



where we assume that the Euclidean space-time rotational symmetry is not retained on the 
lattice. The coefficients Z A U and c^' ' H ' L ^ are functions of the quark masses rriQ and m q . 
With the aid of equation of motion we are allowed to set c Ao = c l Aq = 0. In the special case of 
mq = m q , c A = and c A = c\ is derived from the charge conjugation symmetry. We also 
note that all the improvement coefficients except c A vanishes in the limit of tuq = m q = 0. 
We determine and c^'~'^' L ^ at the one-loop level for both heavy-heavy and heavy-light 
cases. 



A. Determination of the improvement coefficients for the axial vector currents 



The general form of the off-shell vertex functions at the one-loop level on the lattice are 
given by 

Agfog.mpi.rnpa) = 7fc75*f + 7*75 + jUf} + {^f+^b^s 



+(Pk ~ Qk) IsGf + 75 ^Gf + i^Gf + ^G 



fc5 



ik5 



fc5 



rfc5 



rk5 



+ 0(a 2 



7fc75< 5 



(59) 



A$(p, q, m pl , m p2 ) = 7o75^i° 5 + 7o75^ 2 ° 5 + ^7o75^3° 5 + ilol^Ff 1 

+(Po - qo) \l,Gf + ^JGf + <t l5 Gf + MG 



+ (Po + qo) [75^i° 5 + 75^ 2 ° 5 + ^75^ 3 05 + MK] + 0(a 2 ), (60) 

where the coefficients F^ 5 , G^ 5 , are functions of p 2 , q 2 , p ■ q, m p \ and m P 2- 

Sandwiching (f5*H|) and (joTlj) by the on-shell quark states u(p) and u(q), which satisfy 
]/>u(p) = im p iu(p) and u(q)tf = im p2 u(q)^ the matrix elements are reduced to 

W(9) A A5 (P: Q, m pl> m p2 )u(p) 
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and 



u(lhkl5u(p) { F i 5 + im pl F^ 5 + im p2 Ff - m pl m p2 F^} 

+u{q)lkl^ s u{p) [Ff + %m p2 Ff + im pl F* 5 } 

+u(q)4s r Ykl5u(p) {F 5 fc5 + im pl F^ + im p2 F^) 

+ (pk - QkMqhMp) {Gf + im pl Gf + im p2 Gf - m pX m p2 Gf } 

+ (pk + QkMqhMp) { H i 5 + irn pX Hf + im p2 Hf - m pl m p2 Hf} + 0(a 2 ), 



u(?)ASjf (p, q, m pU m p2 )u(p) 

u(qholMp) { F ? 5 + im p iF° 5 + im p2 F^ - m pl m p2 F° 5 } 
+ (Po ~ q )u{q)u(p) {Gf + im pl Gf + im p2 Gf - m pl m p2 Gf} 
+ {Po + qoMq)lMp) {Hi 5 + irripiHf + im p2 Hf - m pl m p2 Hf) + Q{ 



a 2 ) 



where the coefficients are summarized as 

X k5 = Ff + im pl Ff + im p2 Ff - m pl m p2 F* 5 , 
Y k5 = Gf + im pl Gf + im p2 Gf -m pX m p2 Gf, 
Z kb = iff + %m pl Rf + im p2 Rf - m pl m p2 Hf, 

Rkt = -F 3 fc5 + im p2 Fj 5 + irripiFg 5 , 
S k5 = F* 5 + im pl F* 5 + im p2 F$ 



and 



^05 = F° 5 + im pl F° 5 + im p2 F° 5 - m pl m p2 F 4 05 , 
Y 05 = Gf + im p iGf + im p2 Gf - m pl m p2 Gf, 
Z 05 = iff + im pl Hf + im p2 Hf - m pl m p2 Hf. 

In terms of these coefficients the improvement coefficients in eg. (|58)l are given by 





— (x k5 


\latt 


- (X fe5 ) cont 


lC A k 


= {Y k5 ) 


latt 


- (n 5 ) cont , 


lC ~A k 


= (z k5 


jlatt 


/ ry \COnt 

— \^kb) 


~ ic A k 


— (S k5 ) 


( latt 
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< 


= (Rk5 


^latt 


i 


A 7 o75 


= (-^05 


^latt 


- (x 05 ) cont 


ic %> 


= (Yob) 


latt 


- (*05) C ° nt , 




= (Zo5, 


latt 


I ry \ COIlt 

— 1^05 J 



(75) 
(76) 
(77) 
(78) 

where we calculate the continuum contributions employing the MS scheme with NDR. It 
should be noted that Rk5 = Sk5 = in the continuum from the space-time rotational 

for — m p2 from the charge conjugation 
symmetry. 

Combining A 7m75 and the wave function renormalization factors, we obtain the renormal- 
ization factor of the axial vector currents: 



symmetry and c L Ah = c Ak and c A 



ylatt 
^cont 



7 (0) / (0)- 



7(0) / (Oh 
■ / g 1 lattV'' t p2 , 



1 - S 2 A , 



with 



Z QMt{ m pi) = cosh(mJ,° 1 ) ) + r t sinh(m, (0) 



(o) / (oh 



cosh(mp^) + r t sinh(m p 2 ; ) 



pi 
(oh 



A, 



(79) 

(80) 
(81) 
(82) 



2 2 

where Aq i9 are found in Ref. |2J. For convenience we give the relation for Z A ont between 
NDR and DRED in MS scheme: 



Z^(NDR) = (DRED) - -<f 



(83) 



The relevant coefficients ^5,^5,^5,^5,^5 for are determined from the off-shell 
vertex function ([59)1: 

1. 



•yd 1 yd 
/c5 ^£5 



r fc5 ^5 



X fc d 5 + 2im vl R d h 



X d k5 + 2im p2 S d k5 



-Tr 



4 



Al 1 5 ) (l + 7 o) 7 57fc 



p=p ,q=q. 



<9 



d 



., A l 1 5 ) ( 1 + 70)75 -^2(1 + 70)7^75 

Ag?(i + 70)75 + ~-a2(i + 70)7^75 



1 j^fc 



9 



Jp=p*,g=gJ 



4 
1 



-Tr 



<9 



-Tr 



A i5 757fc(l - 7o) + 2im p i— Aj^(l + 7o)7i7s7fc 
A^l + 7o)757fc + 2im p2 T^-Ak\l + 7o)757fc7i 



p=p*,q=q* 
i^k 



p=p*,q=q' d 
i^k 

p=p*,q=q* d 



(84) 
(85) 

(86) 

(87) 

(88) 
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where p* = (p = im pl ,pi = 0) and q* = (q = im p2 ,qi = 0) or q* d = (q = -im p2 ,qi = 0). 
We obtain X 05 , Yqs, Z 05 for A from eq. (joT?j) : 



X, 



05 
1. 



-Tr 



AS 757o + im p i— A$(l + 7o)7fc75 + im p2 —A { ^{l + 70)7*75 
CP/t oqk 



lp=p*,q=q* 



X S 5 + im pl (y o fl 5 + Z* 5 ) - im p2 (Y % - Z ( 



05 J 



-Tr 



d 

Ag 75 (l + 70) + 2im pl — AS(1 + 7o)7fc7s 



< - zm pl (F d 5 + Z* 5 ) - im p2 {Y* 5 - Z, 



05 J 



-Tr 



p=p*,q=q* 
(*)/ 



A 05 j (l + 7o)75 



p=p ,q=q d 



,(89) 

(90) 
(91) 



where F , G and if are functions of p , g and p • q resulting in 







QPH5 


dpi 


p=p*,q=q* 


dqi 






dH^ 


dpi 


p=p*,q=q* 


dqi 








dpi 


p=p*,q=q* 


dqi 



0. 



p=p* ,q=q* 



p=p*,q=q* 



p=p • ,q=q 



(92) 
(93) 
(94) 



with i = 1,2, 3. 



As for the counterterm to isolate the infrared divergence in the above coefficients at the 
one-loop level we employ 



9 (A 2 - k 2 )C F i-f, 



1 



(0) 7^75 



p2 



11a 



k 2 + X 2 



(95) 



with a cutoff A (< ir). 



B. Results for the improvement coefficients of the axial vector currents 



1. Heavy-heavy case 

With the choice of m pl = m p2 the improvement coefficients A 7fc75 , c A ' H ^ and A 7075 , c\ o are 
determined from eqs. (|7T|) . (|72^) and (J7oT -l77 |) . It is noted that c L Ak = c Ak and c~ Ak = c Aq = 
from the charge conjugation symmetry. The quark mass dependences of A 7fe75 ,c^'^^ and 
A7075 > c a are shown in Figs. |H1 and respectively, employing the plaquette and the Iwasaki 
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gauge actions. We also give the dependence of in Fig. [TUl The solid lines denote 
the interpolation with the use of the following functions: 



A 



7fc75 



A* 



i 7075 



4 



Aa 



with a n - and a d l 



A 



+ 



m«?=0 



+ 



Ej=l Q nJ{ m pl' > } t 

l + E^a^mg^ 



1 + Ef =1 ai + W} 



E 



r (0)i i 



i + E 

A A 1 



=1 a di 



(°) n + 



,(°) 
1 

E 5 



i=l "Vm I '"pi J 



l + E 



/ 1 {^pi y 



A 



7075 I 



^=0 + 



+ 



Ej=l a ni{ m pl y 

Ef 



-<i=l a ni 



,0-7, 
ft 1 

(c 
'-pi 



0+ {m ( Sy 



Ef=l a nf { m pl^} 1 



A A 



(°) n + 



1+ELiaH^gT 
5) given in Tables IS and |^ 



of these interpolations are a few % over the range < m. 



pi 



(96) 
(97) 
(98) 
(99) 
(100) 
(101) 
(102) 

The relative errors 
< 10. The mass- 



less values are A. 



7fc75 



A 



7075 



=0.03873(plaqette), 0.04184(Iwasaki), 0.04359(DBW2), 



=0.0075 74(plaquette), 0.003801 (Iwasaki), 0.001492(DBW2) and A Ak 



A Ao =0.1165(plaquette), 0.05663(Iwasaki), 0.02330(DBW2). The coefficient c H A should van- 



ish at m 



(0) 

pi 



2. Heavy-light case 

The set of improvement coefficients are determined from eqs. (j71H 78 )) . For the same 
reason before we keep = 0.0001 for the determination of c\ k and c^~\ while the 
vanishing light quark mass = is employed for other coefficients. 

Numerical results for A 7fc75 , c^' ' H ' L ^ and A 7075 , c^' ^ are presented in Figs. HP and fT2l 
respectively, in the case of the plaquette and the Iwasaki gauge actions. We also plot the 
7T7.pi dependence of A^ M in Fig. [T3J The solid lines represent the interpolation expressed as 

A,„ = A T „ L ^„ + (103) 
14 



\ 

- A 7075 

Aa 



El 



+ 



E 



1=1 



a^{mg ) } i 



1 a m 



k(-,H,L) 



l + E 



5 

i=l 



h <-^{mg>}* ! 



A 



7075 I 



(0) 

pi 



+ 



ELi4 z K°iT 

Ei=l a nl{ m pl' 1 }* 



l + E 



°A 

n 



+ 



E 



=1 u nt 



;1 a ni {Wpi'*} 1 



1 + Ef=i 



{-ill 



(0)- 
pl . 

Ef=l a nf {^pl''}' 



l + E 



OZ 



{^ P l} 4 



where the relative errors to the data are a few % over the range < rn$ < 10 



(104) 
(105) 



(106) 
(107) 
(108) 
(109) 



of the parameters ajjf and a^-° (i 



the constraint that c\^ n ^' and c^ vanish at = 0. 



The values 

5) are listed in Table IVHl and IVTTT1 Here we use 

,(°) 



IV. MEAN FIELD IMPROVEMENT 



Let us explain the mean-field improvement on the renormalization factors of the vector 
and the axial vector currents in eqs. ()25|) and (|79|). We first rewritten their expressions as 
follows: 

Z 



^cont 



^latt 

^cont 
A a 



latt 



+S ~2~ MF + 2 Z (0) gm (0) + 2^T^W Amp2 I ' (110) 

Z Q,latt am pl Z g,latt °' m p2 



, 2°F T J- # aZ/ Q,latt a . 1 P OZj qMt A I /n^ 

+s ^-t mf + 2pr-^r A "v + 2 Z w g m (o) Am ^ 2 1 ' (lj 1 ' 

Z Q,latt am pl Z g,latt °' m p2 
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where and Am pl)2 are defined in Ref. |2j and T MF is the one-loop correction to the 

mean-field factor defined by 

u = pV4 = 1 _^CV Tmf (n2) 

with P the plaquette. We find a detailed description on the derivation of Tmf in Sec. Ill of 
Ref. Q. We then replace u by P 1 / 4 measured by Monte Carlo simulation. 

The mean-field improved MS coupling fi'fjg(Ai) at the scale /i is obtained using the lattice 
bare coupling g\ and P: 

IP 22 / 4 \ 

g^y = 7o + da + Cp + logM + Nf r ~ w log W J (113) 

with iVy the number of quark flavor. The values of c p , d g and df for various gauge and 
quark actions are summarized in Ref. For the improved gauge action one may use an 
alternative formula jl^ 

1 c P + SctRl + 16c 2 P2 + 8c 3 P3 

22 

+d g + (c • c p + 8ci • c m + 16c 2 • c R2 + 8c 3 • c m ) + log(/«i) 
+^/(rf/- i ^logM), (114) 

where 

(115) 
(116) 
(117) 
(118) 

and the measured values are employed for P, PI, P2 and P3. We also find the values of 

, n 

crx, cr2 and cm for various gauge actions in Ref. |11| . 



p = 


l^u pl = 1 - 


-c p gl + 0{g%), 


PI = 


^TrP rt5 = 1 


- c Rx gl + 0(g& 


P2 = 


^TrU chr = 1 


- c R2 gl + O(go 


P3 = 


-TrU plg = 1 


- c R3 go + 0(g%] 



V. CONCLUSION 



In this paper we have determined the 0(a) improvement coefficients of the vector and 
the axial vector currents in a mass dependent way at the one-loop level. Our calculation 
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is made employing the relativistic heavy quark action, which we have recently proposed, 
with the various gauge actions. The results are presented both for the heavy-heavy and the 
heavy-light cases as a function of heavy quark mass. 

For convenience we have given a brief description about the implementation of the mean 
field improvement for the renormalization factors. We are now performing a numerical 
simulation for the heavy-heavy and heavy-light meson systems using the relativistic heavy 
quark action with the 0(a) improved vector and axial vector currents, whose parameters 
are mean-field improved at the one-loop level. This work would reveal to what extent our 
relativistic heavy quark formulation is quantitatively efficient to study the heavy quark 
physics. 
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CO 



TABLE I: Values of parameters v^ 1, '^' H 
eas. (|44l - ITfl) . respectively. 



and Vjp' + ' H ' Z ^ 



(i = 1, . . . , 5) in the interpolation of A 7fc , Cy k ' , Ay fc for heavy-heavy case with 



gauge action v, 



7)1 



J n2 



Vn3 



v n5 



V. 



dl 



■d2 



' ; J3 



V 



do 



A 



a- 



plaquette 

Iwasaki 

DBW2 



0.014684 

-0.015146 

0.010395 



0.13219 

1.8778 

-0.047110 



0.89527 

-0.49846 

-4.5545 



-0.11174 

-1.5767 

0.89861 



0.31131 

-0.31618 

-1.0730 



7.3765 
169.33 
63.029 



34.893 
117.23 
44.429 



2.4179 
127.59 
3.2276 



14.916 
13.898 
13.901 



0.070066 

0.0027264 

-0.022777 



'V k 



plaquette 

Iwasaki 

DBW2 



25.076 

0.94447 

-0.0022351 



-181.54 
-8.8867 
-0.00014326 



-195.34 

2.5958 

-0.034665 



2.4509 
23.741 
-0.043163 



-4.1924 
-1.9005 
0.047341 



22843. 
1351.1 
-5.2267 



7808.2 

-2317.0 

20.895 



13344. 
361.82 
-17.228 



-324.79 
-2370.9 
16.897 



260.72 
187.06 
-9.3839 



Ji 

'Vl. 



plaquette 

Iwasaki 

DBW2 



0.037455 

-0.0029594 

-0.088565 



0.14443 

-0.011657 

-0.18583 



0.32618 

0.0030631 

-0.22715 



-0.029729 

-0.25435 

-0.019377 



0.056502 

-0.37900 

-0.19151 



12.703 

-1.3206 

6.4502 



14.177 
9.8642 
6.1322 



41.911 
51.400 
7.0131 



3.7738 
29.745 
2.7602 



6.8538 
31.961 
4.6325 



Ay, 



plaquette 

Iwasaki 

DBW2 



-0.16445 

-0.057093 

-0.026515 



-7.8793 
-10.029 
-1.6401 



0.59083 
-5.0701 
-5.5463 



-0.45059 
-0.97000 
-1.0198 



0.058239 

-0.56750 

-2.0489 



42.738 
116.04 
43.104 



36.668 
153.83 
62.938 



2.0995 
4.2651 
17.548 



3.3521 
24.187 
24.974 



0.0079208 

-0.027538 

-0.013591 



TABLE II: Values of parameters 
eqs.(@HH2I31)) respectively. 



and v 



o(rr,+,z) 

di 



(i = 1, ...,5) in the interpolation of A 7o ,Cy , Ay for heavy-heavy case with 



to 
o 



gauge action v. 



o 



J n2 



V 





;i3 



„0 
v nA 



,0 

-«5 



J dl 



,0 

'(12 



V 





</3 







,.0 

'c/5 



plaquette 
A 70 Iwasaki 
DBW2 



0.00019720 2.7579 -4.4172 -0.95205 -0.047970 90.897 50.506 33.105 0.57994 0.021081 
0.015673 -1.1138 5.9122 -6.4339 -3.4423 -64.281 99.711 73.380 66.707 -0.35366 
0.012159 -0.021725 -4.5857 -0.90786 -0.89745 38.421 46.673 14.201 9.0663 -0.0099410 



plaquette 
Cy o Iwasaki 
DBW2 



0.0071689 

0.019308 

-0.098707 



-0.34738 0.011308 
-4.4339 -1.9827 
-0.67466 -1.2994 



0.042259 

0.50641 

-0.58528 



-0.17224 
-0.19121 
-0.047428 



15.931 
148.80 
9.5348 



-6.9566 

57.869 

25.681 



12.093 
62.367 
22.818 



-2.5858 
-13.529 
10.537 



4.7195 
4.7863 
0.82699 



plaquette -0.16426 -9.6984 -3.3551 -1.0849 -0.26984 56.744 44.807 5.8125 6.4396 0.0055241 
Avb Iwasaki -0.077730 -1.6754 -6.3923 -1.1045 -0.99597 21.914 79.398 4.6312 16.091 0.0065184 
DBW2 -0.026338 -0.76017 -3.1685 -1.5062 -0.78087 18.143 29.877 13.287 7.5750 -0.0039958 



TABLE III: Values of parameters v. 
case with eas.(|51l-(M|). respectively. 



k(rf,+,-,H,L,Z) 



and v 



k{f,+,H,L,Z) 
di 



(i = 1, . . . , 5) in the interpolation of A 7fe , c 



(+-,H,L) 



Ik > ^Vk 



, Ay fe for heavy-light 



gauge action v. 



nl 



V. 



n2 



u n3 



«5 



d2 



d?> 



J dl 



plaquette 

Iwasaki 

DBW2 



0.0062511 
0.0031637 
0.0021629 



0.052107 

0.26026 

0.025944 



0.55038 
0.55065 
-0.039416 



0.12282 

-0.20725 

0.14083 



0.081216 
0.060822 
0.045331 



8.3279 
64.787 
4.9552 



94.574 
116.19 
-3.1981 



36.477 

-8.9718 

26.268 



15.865 

0.47615 

25.127 



2.8063 
3.3216 
4.2530 



plaquette 

Iwasaki 

DBW2 



0.0015687 

0.70482 

-0.0053483 



-0.011647 

-8.0596 

0.030398 



-0.56066 

-9.4795 

-0.17270 



0.12164 

-0.071982 

0.20455 



-0.21315 

0.023359 

0.077260 



14.953 
1422.3 
-4.8436 



39.726 
920.72 
23.308 



29.119 
921.97 
-21.287 



5.8968 
5.9010 
-12.747 



12.722 

-2.2079 

-14.273 



plaquette 

Iwasaki 

DBW2 



-0.0022183 

0.0072349 

0.062422 



0.091343 
0.033385 
0.40457 



-0.81318 

0.090949 

-0.18140 



1.0896 

0.14881 

0.15221 



0.025931 

0.0022120 

0.0031810 



-23.899 

3.5137 

15.041 



120.15 
32.527 
49.001 



188.77 
29.132 
-2.6721 



242.08 
42.756 
11.963 



75.209 
8.4052 
7.7052 



Ji 



plaquette 

Iwasaki 

DBW2 



0.0054262 

-0.00071764 

0.025927 



-0.11193 

0.00044474 

-0.65927 



0.67635 

0.020620 

2.7483 



-1.0049 

0.0011663 

-9.1542 



-0.028206 

-0.000070314 

-0.11381 



-25.354 
-1.3698 
-371.59 



249.89 
16.036 
1369.9 



-910.86 

3.7954 

-3734.8 



-443.70 

3.5247 

-1534.4 



-147.73 

-0.088466 

-1256.5 



-V k 



plaquette 

Iwasaki 

DBW2 



0.039461 

0.0035068 

-0.10580 



0.16335 

0.013396 

-2.3541 



0.17557 

-0.0074426 

-3.0241 



0.53702 

0.092313 

-0.75090 



0.31432 
0.10210 
-0.28242 



7.9009 
3.5571 
41.820 



6.4773 

-2.5954 

151.64 



24.665 
24.803 
105.43 



19.259 
16.574 
34.058 



3.7410 
3.6739 
9.0980 



A 



plaquette 

Iwasaki 

DBW2 



-0.092742 
-0.043329 
0.021307 



0.18493 

2.6134 

2.2956 



-1.7660 
-2.8159 
0.36278 



0.073087 

0.29050 

-1.9582 



-0.044986 
-0.061514 
-0.26512 



-1.2430 
-59.887 
-171.10 



16.954 
18.652 
10.537 



13.695 
40.112 
37.125 



1.3841 

-0.18607 

83.567 



0.61047 

1.0233 

2.1427 



TABLE IV: Values of parameters v, 
eas- H54l - IS7j) . respectively. 



0(7r 
ni 



and v 



0(-y,+,-,Z) 
di 



{i = 1, . . . , 5) in the interpolation of A 



c (+ ' _) 

70 > °V ' 



Ay for heavy-light case with 



to 
to 



gauge action 
plaquette 

A 7o Iwasaki 
DBW2 
plaquette 

Cy o Iwasaki 
DBW2 
plaquette 

Cy Iwasaki 
DBW2 
plaquette 

Ay Iwasaki 



v nl 



v° n2 



v n3 



,,0 
dl 



V 




(12 



</3 



V 









</5 



0.012802 

0.0077753 

0.0063022 



0.49100 

0.061860 

-0.083313 



0.33452 

0.024810 

0.56176 



0.23128 
0.14301 
-0.17488 



0.040371 
0.069308 
0.017591 



38.224 
6.9955 
-11.399 



37.290 
7.3898 
87.354 



48.752 
13.856 
7.4716 



16.627 
22.055 
-8.3999 



1.3978 
2.9415 
1.1371 



0.0048159 

-0.0094956 

-0.047989 



-0.022115 

0.12375 

-0.42990 



-0.65883 

0.53923 

-0.37339 



0.0040211 

0.36455 

-0.55380 



-0.17926 

0.29355 

-0.33072 



23.777 

-14.207 

13.176 



68.425 

-47.045 

23.031 



36.502 

-49.638 

24.888 



19.322 

-38.991 

25.931 



9.7039 

-14.935 

10.361 



-0.014388 

0.015770 

0.12368 



0.00032644 
-0.085430 
1.3600 



0.13076 
0.55703 
-0.071861 



0.046446 

0.30343 

0.088198 



0.0015886 

0.022609 

0.14040 



7.5778 

-4.5521 

24.213 



21.438 
32.649 
35.781 



13.630 
32.280 
4.4921 



9.7527 
25.481 
1.9685 



0.74691 

2.2470 

5.1811 



DBW2 



-0.092162 

0.18102 

0.051388 



0.56122 

44.773 

4.0357 



-0.53593 

-167.94 

2.3373 



0.082763 

16.819 

-1.5950 



-0.054690 

-4.2353 

-0.26105 



-5.1290 
-1236.9 
-335.95 



-0.12544 

3747.4 

-109.32 



5.3433 
2472.3 
-56.028 



-1.1662 

43.654 

89.233 



0.88503 

76.008 

2.1341 



to 

CO 



TABLE V: Values of parameters CL k 7 ^ l,+ ' H,Z ^ and o k J[ <,+ ' H ' Z ^ {i = 5) in the interpolation of A 7fc75 ,c^~ . 

with eas.(|§6T- l§§)) . respectively. 



AA k for heavy-heavy case 



gauge action 



a nl 



a n2 



a n3 



a n4 



a n5 







dl 



a 



(12 



l d?> 



a 



dl 



l diS 



plaquette 
A 7fe75 Iwasaki 
DBW2 



-0.019142 
-0.011685 
-0.0045598 



-0.047713 

0.10626 

0.11534 



-0.0062276 

-0.35238 

0.022512 



-0.017168 

0.26829 

0.11025 



0.035701 

0.20572 

-0.00049818 



1.7330 

-7.0320 

5.1886 



5.6259 
16.561 
0.97387 



1.2548 
1.8063 
4.8579 



4.7507 
22.443 
0.79928 



0.27302 

1.3801 

-0.012343 



plaquette 



c"t Iwasaki 
DBW2 



0.018067 

-0.029644 

-0.15917 



-0.14688 

0.010150 

-1.5545 



-0.12254 
-0.27996 
-2.5459 



-0.38779 
-0.30399 
-1.1665 



-0.20974 
-0.12102 
-0.045035 



13.975 

0.51603 

18.164 



9.5264 
9.4466 
67.405 



35.820 
18.071 
67.406 



13.366 
11.317 
33.382 



11.765 
5.4941 
1.2761 



plaquette 
cH Iwasaki 
DBW2 



-0.021064 
0.00045423 
0.072677 



0.11720 
0.21074 
0.62160 



-0.0045585 

0.23023 

0.53520 



0.037022 

0.10905 

0.19668 



0.025896 
0.061551 
0.019528 



7.7376 
15.102 
17.971 



1.2405 
14.308 
38.152 



8.1200 
23.396 
24.444 



0.0044767 

4.7100 

10.407 



2.6000 
5.0676 
1.0515 



plaquette 
A,4 Iwasaki 



DBW2 



-0.17057 
-0.10682 
-0.040023 



-8.9210 
-2.6283 
-0.033273 



-5.1962 
-4.7827 
-0.28351 



-0.59753 
-0.20855 
-0.0094526 



-0.30081 
-0.23399 
-0.010229 



39.166 
26.635 
3.7131 



51.312 
58.249 
6.7739 



35.011 
71.704 
20.057 



1.9229 

0.47394 

-1.3448 



2.9460 
5.7885 
1.6373 



TABLE VI: Values of parameters a 
eas. ([T0Ul - ITU21t . respectively. 



ni 



and a^[ 1,+,Z ^ > (i = 1,...,5) in the interpolation of A~ /0 ,c\ o , Aa for heavy-heavy case with 



o 



gauge action a„ 



a n2 



7° 
'ii3 



7° 



a 



o 

71 5 



7° 

VI 



7° 

V2 



7° 

V3 



7° 



a 



o 

(15 



plaquette 0.0012408 -0.027454 0.15598 
-0.00032900 0.041856 0.49428 
DBW2 0.0024558 0.10332 2.0082 



A 7075 Iwasaki 



-0.17366 0.55058 -7.9684 
-0.078257 0.24277 12.589 
-0.81526 0.40598 13.920 



80.466 -14.079 41.179 3.5040 
29.569 -4.1770 12.693 1.4849 
59.897 -15.834 9.6390 2.3878 



-A 



plaquette 0.052844 0.098447 -0.095352 0.000049485 0.00013520 13.949 
Iwasaki -0.020399 -0.22848 0.42451 -1.5113 -0.93644 14.179 
DBW2 -0.21817 -4.0728 1.8993 -2.2319 -0.061938 33.245 



9.2995 7.5620 4.0175 -0.050684 
-31.091 111.08 36.737 58.550 
78.786 -31.755 49.858 1.5508 



plaquette 
A^ Iwasaki 
DBW2 



-0.19696 -8.7519 -5.3415 -0.48001 -0.10054 44.563 
-0.094508 -2.4786 0.35337 -0.032386 -0.00029766 28.057 
-0.035180 0.032360 -0.029443 0.0065747 -0.00021265 2.8063 



55.001 31.792 5.2293 1.3789 
21.317 4.6126 -0.66769 0.11312 
-2.8195 3.0682 -0.27495 0.033404 



AC- 



TABLE VII: Values of parameters a, 
case with eas. (jl03l - ll06j) . respectively. 



,H,L,Z) 



and a 



k(y,+,H,L,Z) 
di 



, 5) in the interpolation of A 7fc75 , c ~^' ,H ' L ^ f A^ fe for heavy-light 



gauge action a 



<n2 



x n3 



l n4 



''7(5 



dl 



l d2 



'd?> 



'■di 



'(15 



A 



7fc75 



plaquette 

Iwasaki 

DBW2 



-0.0077994 

0.010623 

-0.0027159 



-0.41731 

-8.7336 

-0.023333 



1.1354 

-0.14090 

0.11288 



0.094127 
-1.3650 
-0.091528 



0.29482 

-0.040416 

-0.0081287 



45.381 
1596.5 
18.209 



-95.269 

1554.8 

-16.459 



-62.935 

1388.6 

-67.226 



-98.936 

375.05 

-42.775 



-29.162 

9.6218 

-14.363 



plaquette 

Iwasaki 

DBW2 



0.0019484 

-0.0085895 

-0.061215 



-0.12572 

0.10146 

0.080878 



0.15951 

-0.15900 

-0.031488 



0.073993 

0.32234 

0.019396 



0.052235 
0.081455 
0.0023197 



-22.185 
-9.2200 
6.6362 



-50.850 
-11.177 
-7.8813 



-46.419 
-1.5297 
1.1771 



-26.709 
-89.575 
-1.5133 



-6.6887 
-19.944 
-1.2521 



plaquette 

Iwasaki 

DBW2 



0.0079671 

-0.0015928 

0.0036783 



-0.37407 

0.010270 

0.00027955 



2.8843 

0.057064 

0.016457 



-5.3964 

0.0053304 

0.0096355 



-0.14091 

-0.00012449 

0.000079022 



27.423 
1.8033 
-1.0696 



238.43 
15.634 
5.2927 



-806.02 

4.7995 

1.5712 



-384.25 

3.5412 

2.3740 



-198.31 

0.0063867 

0.31191 



Ji 



plaquette 

Iwasaki 

DBW2 



0.0017201 
0.0010232 
0.0010612 



-0.0031462 

-0.000089135 

0.0021013 



0.053022 
0.030173 
0.014049 



0.0013268 

0.00049781 

-0.00012538 



-0.000064448 
-0.000051947 
-0.000040511 



-1.1603 
-1.2988 
-1.6824 



14.397 
11.412 
8.7455 



4.2428 
2.6448 
0.27197 



3.3529 
2.7244 
2.1129 



-0.040156 
-0.067972 
-0.11394 



plaquette 

Iwasaki 

DBW2 



-0.028427 

0.017833 

0.092039 



-0.12517 
-0.64125 
0.66663 



-0.025615 

-3.1710 

-0.51025 



-0.15291 

4.4169 

1.1278 



0.010900 

0.64578 

0.22717 



10.057 

-242.31 

16.138 



23.991 

-174.96 

5.1945 



10.071 
214.38 
9.1509 



34.909 
450.52 
32.237 



3.2330 
33.050 
4.3118 



plaquette 

Iwasaki 

DBW2 



-0.10998 

0.12336 

-0.11610 



0.30354 

-42.429 

-12.061 



-2.9348 

3.5897 

-3.5918 



0.12711 

-2.2117 

-1.7161 



-0.084194 
-0.054401 
0.0095736 



-1.8757 

773.10 

621.73 



23.675 
590.97 
304.15 



19.521 

-12.310 

192.74 



1.4942 
42.384 
23.838 



0.75518 
0.13072 
-0.43073 



to 



TABLE VIII: Values of parameters a 
with eas. (|lU6l - llU9j) . respectively. 



0(7r 
ni 



,Z) 



and a°- 7 ' + ' ' 



(i = 1, . . . , 5) in the interpolation of A 



,(+-) 



7075 ' Aq 



, Aa for heavy-light case 



gauge action 



a nl 



7° 



*n3 



a n4 



a n5 



a 







7° 

V/2 



7° 



a 



o 



a 



o 

</5 



A 



7075 



Aa 



plaquette 

Iwasaki 

DBW2 



0.000068956 
-0.0032698 
0.0016117 



0.00047409 

0.046507 

-0.018932 



0.037244 
0.085071 
0.27891 



-0.015986 
-0.098217 
-0.16669 



0.013863 

-0.0091805 

0.060548 



4.6946 
54.923 
20.906 



21.769 

-14.421 

58.558 



-6.2023 
-17.932 
-19.071 



6.2455 

-10.622 

3.4284 



0.67199 

-0.41183 

4.5131 



plaquette 

Iwasaki 

DBW2 



0.016913 

-0.012639 

-0.10453 



0.098293 

0.12162 

-0.42038 



0.13800 

-0.20061 

0.25130 



0.11672 

-0.49574 

-0.44821 



0.058833 
0.094685 
-0.11844 



9.2871 

-8.5494 

10.887 



49.110 
5.8278 
7.7564 



25.152 
52.240 
2.5891 



12.626 
44.235 
13.146 



1.7409 
10.945 
6.4704 



plaquette 

Iwasaki 

DBW2 



-0.011811 
0.0071466 
0.046117 



-0.035056 

-0.24087 

0.28507 



-0.57156 

0.15744 

-0.026074 



-0.17770 
-0.21666 
0.24895 



-0.36685 

0.13313 

0.045710 



5.9258 

-47.112 

10.757 



76.946 

-65.362 

10.777 



38.435 

-38.220 

6.7615 



42.946 

-50.492 

11.691 



8.8915 

-10.153 

2.7523 



plaquette 

Iwasaki 

DBW2 



-0.10088 

-0.024699 

-0.011150 



0.73718 

-6.5964 

1.7357 



-4.0136 

-0.47808 

-1.4097 



0.12455 

-0.45893 

-0.44819 



-0.11617 

-0.018953 

-0.13594 



-6.4483 

136.12 

-98.330 



33.715 
130.02 
54.859 



29.058 
14.889 
32.635 



2.8742 
13.667 
28.748 



1.3900 

0.091530 

1.5363 




FIG. 1: One-loop diagrams for the vertex functions, q denotes the outgoing quark momentum and 
p denotes the incoming quark momentum. 
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FIG. 2: A 7fc , for heavy-heavy case as a function of . Solid symbols denote the plaquette 

gauge action and open ones for the Iwasaki gauge action. 
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FIG. 7: Ay M for heavy-light case as a function of . Solid symbols denote the plaquette gauge 
action and open ones for the Iwasaki gauge action. 
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FIG. 11: A 7fc75 ,c^' ' H ' L) for heavy-light case as a function of . Solid symbols denote the 
plaquette gauge action and open ones for the Iwasaki gauge action. 
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FIG. 13: for heavy-light case as a function of m^. Solid symbols denote the plaquette gauge 
action and open ones for the Iwasaki gauge action. 
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